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The zero gap surface states of a 3D-topological insulator host highly mobile Dirac fermions with 
spin locked to the momentum. The high mobility attributed to absence of back scattering is reduced 
in presence of impurities on the surface. In particular, we discuss and compare scattering times for 
localized impurities on the surface, scattering between states of opposite helicity located on different 
surfaces coupled through a hybridization potential, and the role of magnetic impurities. Magnetic 
impurities give rise to an additional spin suppression factor. The role of warped bands and its 
influence on topological factors that can enhance the overall surface mobility is examined. 


I. Introduction 

The zero gap helical surface stateJi^ of a 3D topological 
insulator (TI) film can be split by an external magnetic 
field or a reduction of film thickness.l^ Such a band gap 
open top olog ical insulator loses protection against back 
scatterinj^ and has a reduced surface mobility. The 
suppression of mobility is also observed through large- 
angle scattering in presence of impurities on the surface 
of a TI with a consequent degradation of its charge trans¬ 
port properties. Further, a sufficiently thin film, with 
preserved time reversal symmetry, allows the possibility 
of a surface electron(a Dirac fermion) to back scatter^ 
by tunneling through the bulk and eventually occupy a 
state of opposite helicity on the bottom surface thus pro¬ 
viding a pathway for velocity reduction. Inter and intra¬ 
band scattering events on the Dirac hyperbola^ of each 
surface, created by splitting the TI states with an exter¬ 
nal magnetic field or a layered ferromagnet offers a way 
to probe phenomena such as weak anti-localization, 133 
electric conductivity, and localized magnetic moments 
through magnetoresistance measurements. 

In this work, within a Boltzmann approximation, the 
impurity scattering time for surface carriers is evaluated 
(Sec. II) in a topological insulator him. The impurities 
are assumed to be non-interacting. Scattering times are 
evaluated for four cases l)a zero-gap non-trivial topo¬ 
logical insulator 2) in presence of inter-surface coupling 

3) with substrate induced intrinsic dipole moments and 

4) with a hnite magnetic moment with an out-of-plane 
component that introduces a spin induced suppression 
factor. The corresponding semi-classical mobilities (Sec. 
Ill) derived from a linearized Boltzmann equation and a 
relaxation time approximation are discussed. The role of 
warpinj^ and its inhuence on overall transport behaviour 
is also highlighted. Results are collected in Sec. IV fol¬ 
lowed by concluding remarks. 

II. Scattering rates for a localized impurity 

Surface states in topological insulators such as Bi 2 Te 3 , 
Bi 2 Se 3 , and Sb 2 Te 3 are characterized by a linear disper¬ 
sion for small values of the momentum vector and a sin¬ 
gle Dirac cone at the T point. The calculations carried 
out in this work are partly performed using a low-energy 


continuum four-band k.p Hamiltonian for 3D TIs. The 
four-band k.p HamiltoniarfiS! is further simplihed in to a 
two-band Dirac Hamiltonian that represents only the lin¬ 
early dispersing surface states; additional modihcations 
that arise on account of higher warping terms are in¬ 
cluded later for comparison. The four-band k.p Hamil¬ 
tonian in the basis set of the four lowest low-lying states 
|P1+ t), \P2- t), \Plt i), and \P2- i) is 


H{k) = €{k) + 


/M{k) 
I Aik^ 

\A2k^ 


Aikz 0 A2k- \ 
-M(fc) A2k- 0 I 
A2k+ M{k) —Aikz I 

0 -Aikz -M{k)J 


( 1 ) 


where e(fc) = C + Dik^ -I- D 2 k\, M{k) = Mq -I- Bik^ + 
B 2 k\ and k± = ± iky. For Bi 2 Te 3 and Bi 2 Se 3 , the 

relevant parameters are taken from Ref. [TTJ The Hamil¬ 
tonian is diagonalized numerically to obtain eigen values. 
The simplified two-dimensional Dirac Hamiltonian for 
the surface states can be written aJ^ 


Psurf .states kvgi^dxky tfykx^ 


( 2 ) 


where Vf denotes Fermi-velocity and cri]i = x,y are the 
Pauli matrices. The corresponding density of states in 
standard notation is 




dk 


{2tt) 

Changing to polar coordinates, 

1 


6{e-E{k)) 


(3a) 


^27r 


D(e) = 


( 27 r)" Jo 

2 

2tt {hvff 


d0 


k dk 


Hvf 


S k- 


hvt 


(3b) 


The density of states of electrons described by the two- 
band Hamiltonian (Eq. with spin degeneracy included 
is linear in energy and vanishes at the Dirac point. 
In deriving this expression, we have used the property 
(5 (—x) = 6 (x). 

This form of the Hamiltonian, though in principle suf¬ 
ficient to probe the surface states, fails to account for 





2 


the underlying crystal symmetries.^1^ A more accurate 
representation (Eq. ^ within the two-band modeP^ that 
is consistent with the point-group symmetry of the 
rhombohedral crystal, the snow-flake Fermi contour ob¬ 
tained from ARPES, and also preserves time reversal 
symmetry must contain higher order terms cubic in k. 
The next set of corrections involve termJi^ which are 
ignored in all calculations here. 

H{k) = €o{k)+Hvf{aj;ky-aykx)+^h^ {k^ + k^) (4) 

eo(fc) introduces the particle-hole anisotropy and the cu¬ 
bic terms denote warping. Using Eq. and ignoring 
particle-hole anisotropy without loss of generality, the 
surface state spectrum is 

£±(^) = hv'jk? + \^h^k^cos'^{W) (5) 

where 6 = tan~^ (ky/kx)- 

The scattering rate on surface of a topological insulator 
due to localized and static non-magnetic impurity with 
a spherical potential V = VqS {r — Rj) can be evaluated 
using the Fermi-Golden rule. The concentration of im¬ 
purities is assumed to be dilute such that there are no 
interference effects between successive scattering events. 
The matrix element is 


where Xkk' — ykj s'j. The additional spin¬ 

scattering factor C ^s, s ^ takes in to account the helical 
spin structure of the TI surface states. For a pristine TI, 
C (^s,s 'j = with 0 ^ ^ TT, since back scattering 

is forbidden and the factor must assume a value between 
zero and unity. The integral in Eq. therefore must be 
integrated over all values of (j). 


- = ^5 (s) [ Vocos'^t (1 - cos(l)) d(t) 

T n Jq z 

^ TTg Vq 

n^vj 4 


( 9 ) 


where we have used Eq. 3b for density of states g (e). 
The impurities on the surface can also be of several types, 
each with a different density, which means that Eq. 
must be modified to reflect this situation 


n, ”'0 


( 10 ) 


where ni and Vbi denote a particular impurity density 
and the related scattering potential. 

B. Scattering rate between opposite surfaces 


( 6 a) 

where the two-component wave function 4* for the linear 
Hamiltonian (Eq. is 


The two opposite surfaces in a topological insulator 
host states of opposite helicity; therefore, for the possi¬ 
bility of scattering to happen between two such states, 
the topological insulator film must be sufficiently thin. 
A thin film TI can be modeled as 
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A. Scattering rate between states of identical 
helicity 

The helicity of the electron wave function is denoted 
by 77 = ± and A represents a band gap for the surface 
states. Evaluating the square of the matrix element for 
a ungapped TI film when an electron elastically scatters 
through an angle (j) between states with identical helicity 
gives 


\M (k',k) p = |(4'/|Uo|4'.)P = Vo^cos^^ (7) 

and (j) = 9f — Oi- The elastic scattering time using a 
Boltzmann approximatiorP^E^ such that |A:| = \k \ is 

1 = yJ ~ ) lAfefe' I (1 - cos(/>) (8) 


^surf.states = hvf{axky - (Tykx) + Ahax (11) 


where 2A;j is the band gap shift. The scattering rate can 
be similarly computed by employing the Fermi-golden 
rule, the square of the matrix element now is given as 


V^siri^ — . Inserting the expression for matrix element in 
Eq. [ 8 j we get 


1 
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^ Uq^ {hvfkf 
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TTE Yq (hvfk)^ 
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Topological insulator thin Aims are usually grown on sub¬ 
strates that render the structure asymmetric by adding 
an intrinsic dipole moment and breaking structural inver¬ 
sion symmetry (SIA). Such a thin film with an inherent 
SIA potential U can be described by the Hamiltonian 


ddsurf,states = hvf{axky - Gykx) + Ahaz + Uax (13) 

The wave functions for the Hamiltonian in Eq. are 


_ 1 f Xy{k)exp{i9) 


(14a) 
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where A,, is 


At, = 1 ± 


^ \J + if>-vfky + Uf + {hvfk^f^ 


(14b) 


and 9 = tan ^ j -— ^ . The corresponding eigen values 

\Kx U) 
are 

e (k) = (^Al + (hvfky + uf + (hvfk^f^ (15) 

The scattering time by a direct application of Fermi 
golden rule in the case of a carrier that tunnels through 
the band gap open asymmetric thin film, for instance, be¬ 
ginning from a +k state on the top surface and settling 
in to a —k state on the lower surface(opposite helicity) is 


^ = ^g{e)VQ^ sin^^cos^ 2 ~ 


en 

Atfivj 


vh 


(16) 


where 7 = 


(^{hvfky + uf' -b (hvfkxY 

(fhvfky + uf + (hvfkxY 


A2. 


C. Scattering due to a magnetic impurity 

Impurities impregnated on the T1 surface which pos¬ 
sess a finite magnetic moment break time reversal sym¬ 
metry (TRS). A broken TRS gives a non-zero expectation 
value for spin polarization along an axis aligned to the 
outward normal to the T1 surface. To compute the elas¬ 
tic scattering rate in this case, we must re-evaluate the 

spin suppression factor Y (s,s^ which was set to cos^^. 
The spin polarization vectors, bearing in mind the helical 
structure of the surface states are therefore given by 


= J) 
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The final spin polarization vector is 

S = Sxi + Syj -b Szk (18) 

The spin suppression factor between two states |fc (0i), s) 
and \k' {62) ,s') suffering an angular elastic ^|A:| = 
scattering of </> = 6*1 — 02 is therefore 
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The magnetic field induced band gap is A = -g^EBz- 

As before, the scattering time can be computed using 
Eq. further, in presence of a magnetic field, time re¬ 
versal symmetry is lost and back scattering is no longer 
forbidden. Carriers on the surface after scattering can 
therefore occupy states of opposite helicity or in other 
words with a reversed k vector (for example, -bfci to —^ 2 , 
such that I fell = I — ^ 2 ! for an elastic scattering). The 
scattering times when the final scattered states have iden¬ 
tical {+ki to +k 2 and |fci| = \k 2 \) or opposite helicity will 
be different and can be calculated as follows. We first 
establish the scattering rate between states of identical 
helicity. 

^ = ‘^g{e) J d(/)|(T+|VbC (s,s ) l^'+)P (1 - cos((i) 

( 20 a) 


Inserting the spin scattering factor using Eq. |19a| and 
wave functions from Eq. | 6 b| and carrying out the straight 
forward algebra yields 
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where k = 


A2 


The equation can be numer- 


A 2 -b {hvfk) 

ically evaluated as shown in Sec III. A similar integral 
for scattering between states of opposite helicity can be 
written 
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III. Influence of warping on surface conductivity 


and 


Under a weak external force, the deviation of the elec¬ 
tron distribution from the thermal equilibrium value can 
be assumed to be small which allows us to linearize Boltz¬ 
mann distribution within the relaxation time approxima¬ 
tion. The conductivity can therefore be written as 




2^,2 


( 21 ) 


e V 


which at T = 0 yields cr = ^ ^ g {E) t, where we substi¬ 

tute for the diffusion constant D in the original Einstein 

vir 

relation a = e^g {E) D as D = for a two-dimensional 

system. The scattering time for each scenario consid¬ 
ered can be inserted from their respective expressions 
derived above. Finally, we note that the expressions 
derived for scattering time and conductivity undergo a 
modification if warping terms are explicitly included in 
the analysis. More precisely, the band gap A modihes to 
A + h^Xk^cosW and the density of states also undergoes 
a change as shown below 
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We have used the identity S (g(a;)) = — Jr- such 


S(x — Xj) 

j \9jxj)\ 


that g{xj) — 0 and there are no multiple zeros. Xj is a 
simple zero of g{x). The function g{x) takes the form 

g(x) = e — ^Jhv'jk'^ + \^ffik^cos'^{W) (23) 

The band gap split A is ignored since it is swamped 
by the warping correction at points in momentum space 


away from the Dirac cones at T. The kj in Eq. 22c are ob¬ 
tained by determining roots of Eq. Expanding Eq. 
yields 


hvjk^ -I- \^frk°cos^(?>9) — = 0 


(24a) 


Rearranging, Eq. |24a| is written as a cubic equation with 
k"^ as the variable 

X^h^cos^{39) (kJ)^ -\- hvjk'^ — = 0 (24b) 

The real solution for is of the form 

(wi-kW2)^^^ - (|wi - a;2|)^/^ (25) 

where loi and a ;2 are defined as follows 

OJi = \( —7777^^ (26a) 
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hvf 
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.V 


27 \XhJcos{39) J 4 \XhJcos{39) J 
Finally evaluating the derivative g'{x), we have 
hv'jk + iX^ff’k^cos^ (39) 


9 '(x) = T 


-I- X'^hPk^cos'^ (39) 


(26b) 


(27) 


Putting all of them together, the density of states at any 
given zero of g(x), xj can be written as 


D(e) = 


I (kvf)^ kJ + X^ffik^cos'^ (39) 

(27r)^ Jo (kvf)^ k -I- 3X^tfik^cos^ (39) 

(28a) 


Simplifying, we get 
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r27r 


d9 k- 


D(e) = 

( 27 r)^ Jo (kvf)^ k + 3X'^hX’k^cos^ (39) 

(28b) 

Noting that energy e is a constant, the final form for D(e) 
is 


r2Tr 


D(e) = 


d9- 


(27)"^ Jo + 2X‘^h^k^cos‘^ (39) 


(28c) 


In simplifying the above expression, we have used the 
Dirac cone energy expression Ehv^k'^ -I- X^ff'k^cos^(39). 
Inserting for k from Eq. |26a[ a nd |26b[ the final ex¬ 
pression for density of states (Eq. 28c) is evaluated nu¬ 
merically in Section IV. 

IV. Results 

All results derived in this paper use the two-band 
model (Eq. [2 ) which employs phenomenological parame¬ 
ters to describe band gap splitting and inter-surface cou¬ 
pling in case of an asymmetric thin film. These param¬ 
eters can be di rectly established from a first-principles 
calculatiorfi^Eo] or obtained, as done here, by numerically 
diagonalizing the four-band Hamiltonian of Eq. The 
required parameters are read off the dispersion plot. The 
Fermi velocity of the surface states is taken as 5 xlO® 
m/s for all calculations shown. 

A. Non-magnetic impurities 

We first plot the dispersion relationships (Fig. for 
a nine and three quintuple layer Bi 2 Se 3 topological in¬ 
sulator film. Each quintuple layer is about I.O nm. The 
three quintuple layer film has its top and bottom surfaces 
coupled (hybridized surfaces) which induces a finite band 
gap. On the surface of a 9.0 nm thick Bi 2 Se 3 topologi¬ 
cal insulator slab which has zero-gap states, a uniform, 
dilute, and non-interacting spherical impurity scattering 
potential of 10.0 meV is assumed to be present. Elec¬ 
trons scatter from such fixed impurities with no internal 
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excitations elastically. The scattering time is computed 
at an energy equal to the Fermi level which is set to 100.0 
meV. At this energy, the density of states (DOS), using 
Eq. 3a is 2.93 x 10“®meV“^A“^. Inserting the DOS in 
the scattering time equation (Eq. yields a value equal 
to 9.076 X 10“^°seconds. It is important to note that the 
Fermi-level can be raised up to 0.2 eV, beyond which we 
reach the bulk conduction bands of Bi 2 Se 3 . The bulk 
conduction bands of the slab can be seen to begin from 
the 0.2 eV mark in Fig. [l] 


> 

Q. 

CO 

D) 


CD 

CO 

c/) 

< 



Film thickness (nm) 



FIG. 1: Dispersion of a 9.0 nm thick Bi 2 Se 3 topological insu¬ 
lator slab around the Dirac point (Fig. [^) at 0.02 eV. The 
dispersion of the thin film (3.0 nm thick)(Fig. shows two 
Dirac hyperbolas when the surface states hybridize. 


FIG. 2: Band-gap opening as a function of BbSea film thick¬ 
ness. A band gap opens because the wave function from the 
two surfaces penetrate the bulk and hybridize. 



Bi terminated surface 


Te terminated surface 


periodic 


For the case of a thin TI Him, the possibility of an 
electron tunneling through the bulk and reaching a state 
of opposite helicity on the second surface must be ac¬ 
counted for. The variation in band gap of the surface 
states for a thin film with varying thickness is plotted in 
Fig.§ Retaining the parameters from the case of scat¬ 
tering between states of identical helicity and setting the 
band gap A in Eq. to 100.0 meV, the corresponding 
scattering time is determined to be 8.309 x 10“^°seconds. 
The scattering time is similar to the one obtained when 
scattered carriers are constrained only to occupy another 
state of identical helicity on the same surface. It is easy 
to see from Eq. and Eq. that the ratio of scatter- 

(hvfk)^ 

ing times between the two cases is - - - 3 -. The 

A2 -F {hvfkf 

ratio is significant when k is close to the F point or the 
Dirac cone but tends to unity at larger k values since the 
hybridization induced band gap split is generally a small 
number. As pointed out later, at points far away from F 
where the Dirac crossing happens, the density of states 
using Eq. is inadequate. 

We next turn our attention to a TI thin film grown 
on a substrate that renders the structure asymmetric by 
creating a spatially-dependent dipole moment. A simple 
way to mimic the effect of a substrate is to model the 
film with two different surface terminations as depicted 
in Fig. The asymmetric potential gives structural in¬ 
version asymmetry (SIA) which means that the two sur¬ 
faces do not have equi-energetic Dirac cones. As shown in 
Fig. HI a TI film with bismuth and tellurium surface ter¬ 
mination possesses a non-zero dipole momentPlHlS! with 
an oscillating charge pattern in bulk of the film. In 


FIG. 3: A 9.0 nm thick BETea ultra-thin body oriented along 
the x-axis.The two surfaces have Bi and Te termination thus 
making them chemically non-equivalent. 


absence of the SIA potential, the surface states consist 
of spin degenerate conduction and valence band^^ sep¬ 
arated by the hybridization gap A while the substrate 
induced asymmetry leads to a Rashba splitting.!^ The 
actual band gap (Eq. 151 varies with the overall dipole 
moment. For the thin film considered here, U is found 
to be 1.0 eV^ and using identical parameters as before, 
the scattering time is 9.025 x 10“^°seconds. In evaluat¬ 
ing the scattering time, the wave vector is assumed to 
only have a k^ component equal to 0.1 1/A. 

The scattering time in all the three cases is roughly 
the same which suggests that the band gap open¬ 
ing and the dipole moment do not significantly af- 

r ^ 1 • 1 r ihvfk)^ , 

feet the results since the factors - - -?r and 


A2 -F [hvfl^ 


from Eq. 


12 


and 


(j^hvfky + U)^ + {HvfkxY 

A2 -F (^{hvfky -F U)^ -F {hvfkxY 
Eq. ^ at moderate k vectors are numerically close. These 
factors, nonetheless, have a relatively important contri¬ 
bution when the ratio of band splitting energy (A) to the 
Fermi-energy hvfk is not negligible which is possible for 
very small values of the momentum vector. 

All the scattering times discussed here include the spin 
suppression factor whose contribution can be gauged by 
evaluating the scattering rate equation for a zero-gap TI 
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FIG. 4: The spatially-dependent electrostatic potential 
(Fig. 1^) and charge on each atomic node is plotted against 
the X coordinate of the Bi 2 Te 3 thin film. This thin-film has 
Bi and Te termination on the surfaces. An oscillating charge 
pattern (Fig. Ep) is obtained for non-equivalent surface ter¬ 
mination. 


with impurities (Eq. 


91 without the cos^ ^ term. Such a 


calculation gives a scattering time of 2.8 x 10“^°seconds, 
which is a reduction by a factor of 3.24. A shorter scat¬ 
tering time signifies lower mobility and conductivitj^^ 
which can be easily explained by noting that the high 
mobility of surface states is attributed to spin-protected 
forbidden back scattering. The condition to forbid back 
scattering is relaxed by ignoring the spin-suppression fac¬ 


tor of cos^ —. We also note, as shown in Ref. EH the 

scattering time in a gapped topological insulator (trivial 
case) is four-fold higher compared to a zero-gapped topo¬ 
logical insulator (non-trivial case). The factor of 3.24 in 
the present case where the suppression factor controls 
back scattering matches well with the previously derived 
result. In passing we mention here, that the dipole mo¬ 
ment, in principle, at a critical value, can close the band 
gap and re-open as an instance of topological phase tran¬ 
sition.!^. 


B. Magnetic impurities and warped bands 

The helical natur^^SMI Qf surface states of a pris¬ 
tine topological insulator ensures that there is a definite 
in-plane spin polarization. Magnetic impurities on the 
surface give an additional out-of-plane spin component!^ 
thus essentially turning the spin-polarization to a three- 
dimensional vector. Scattering times corresponding to a 
9.0 nm and 3.0 nm thin film are presented for discussion. 
For the 3.0 nm thin film, we consider the case of a carrier 
tunneling to the other surface and occupying a state of 
opposite helicity. Proceeding along same lines as for non¬ 
magnetic impurities with the spin suppression factor now 
modified such that spins are aligned to the effective spin- 
polarization vector, yields scattering times equal to 7.530 
X 10“^°seconds and 6.935 x 10“^°seconds for the 9.0 nm 
and 3.0 nm TI films respectively. As a first check, the 
scattering times are shorter than 9.076 x 10“^°seconds, 
which is the corresponding result for a zero-gap topolog¬ 
ical insulator film. The presence of a band gap reduces 
the scattering time, augments the large angle scattering 


X 10"® 



FIG. 5: The density of states for a warped dispersion devi¬ 
ates from the expression obtained for linear bands at higher 
energies. At higher energies, the warping terms signihcantly 
contribute to the overall energy spectrum of the TI surface 
states. The warping strength in this calculation was set to 
500.0 evA-T 


and consequently lowers the mobility as expected. In 
calculations shown here, we have explicitly neglected the 
orbital coupling of the magnetic field through a Peierls 
substitution and retained only the spin coupling through 
the Zeeman term. The Zeeman term can be attributed 
to an externally applied magnetic field or the exchange 
field of a ferromagnet layered on the surface of a TI.E2 
At this point, it is worthwhile to establish the change 
in scattering rates due to warped nature of bands, if the 
momentum vectors are not in the immediate proximity 
of the r point. The warped bands modify the density 
of states as derived in Eq. |28c[ A comparative plot of 
density of states obtained from a linear dispersion and 
with warping terms included is shown in Fig. The 
density of states for the warped case was determined by 
numerically integrating Eq. |28c| 

The scattering time for a carrier described by a warped 
Hamiltonian (Eq. on surface of a zero-gap topological 
insulator whose initial and final states are of same helic¬ 
ity is 1.541 X 10“®seconds. The scattering time is higher 
compared to 7.530 x 10“^°seconds because of a lower den¬ 
sity of states with warped bands. It is well-known from 
semi-classical transport theorjEH that a lower density of 
states reduces scattering and increases mobility. Scat¬ 
tering times for a range of band gaps for the linear and 
warped dispersion models is shown in Fig. The scat¬ 
tering time decreases almost linearly under higher mag¬ 
netic fields(increased band gaps) for both the Dirac and 
warped dispersion which is indicative of a greater mag¬ 
netoresistance. The origin of the enhanced magnetoresis¬ 
tance is through a large angle scattering (1 — cos9) term 
that enters the analysis in the scattering rate equation. 
Further, since the electric conductivity, as evident from 
Eq. is directly proportional to the scattering time, the 
warped model would yield a higher value. In other words. 
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A (meV) 

FIG. 6: The scattering time on surface of a magnetically 
gapped topological insulator is higher for the warped band 
model (warping strength = 500.0 evA“® ) because of a lower 
density of states. The scattering time also decreases with a 
greater magnetic field induced band gap indicating lower mo¬ 
bility or a higher magnetoresistance. 


electric conductivity is an implicit function of momentum 
space; at k points away from F, the warped Hamiltonian 
which is a better description of electronic spectrum of 
surface states yields lower density of states and a higher 
electron velocity. 

V. Surface conductivity and mean free path 


The conductivity on the impure surface of a topologi¬ 
cal insulator can be estimated by a direct application of 
Eq. |21[ We consider a non-magnetic impurity of density 
rii = 0.5 X 10^°. The scattering potential is 500.0 eV. The 
scattering event takes place at a Fermi level set to 200.0 
meV for which the density of states (using Eq. 3a for 
linear bands) is 5.877 x 10“®meV“^A“^. The scattering 
time for this case turns out to be 9.076 xl0“^^ seconds. 
Putting all these numbers in Eq. the surface conduc¬ 
tivity, expressed in terms of the quantum conductance 
2e2 

a = —— is 137.595cr. The conductivity at large k values 
h 

is reduced since the corresponding density of states must 
now be computed with warped bands. The DOS with 
warped bands at 200.0 meV evaluates to (by numerical 
integration of Eq. 28c) 4.51xlO“®meV“^A“^ while the 
scattering time is 1.183 xl0“^^ seconds. The surface 
conductivity value of 137.721cr is only marginally higher 
than the corresponding 137.595 a for linear bands. An¬ 
other useful that can be easily computed is the mean free 
path given as VfT. The numbers for linear and warped 
bands are 45.38 nm and 59.15 nm respectively. 

While the above analysis does not account for scatter¬ 


ing events which involve phonons and charged impurities, 
the mean free paths obtained suggest that transport in a 
topological insulator channel material for a conventional 
field-effect transistor in a miniaturized device will be bal¬ 
listic. The experimentally determined mean free path in 
GaAs, a high-mobility semiconductor, is approximately 
34.0 nm at room temperature. 

VI. Conclusion 

In this work we have theoretically presented the evalu¬ 
ation of scattering times on the surface of a 3D topologi¬ 
cal insulator using Boltzmann transport. The role of spin 
suppression factor which comes from the topological na¬ 
ture of the states and explains the higher mobility on a TI 
surface is highlighted. It is found that in case of a mag¬ 
netically split TI surface, the alignment of the spin polar¬ 
ization vectors to the effective magnetic field reduces the 
overall scattering time. The influence of density of states 
from a warped model is analyzed and a higher mobility 
and electric conductivity value is established. The en¬ 
hancement is attributed to the reduced density of states. 
In carrying out these calculations we have assumed a uni¬ 
form Fermi velocity for the surface states though that 
may not be true for all materials that exhibit topological 
insulato r beha viour. In particular, the topological Kondo 
insulatorl^^*^ which has Fermi pocket^^with varying ve¬ 
locities is a case in point. Also, in a more careful calcu¬ 
lation, which could be potential future work, preferential 
scattering directions dependent on the underlying crys¬ 
tal symmetry must be considered. Finally, we note that 
because of the helical structure of the surface state, the 
spin polarization vector must have a definite polarization 
for a given k state. It is tacitly assumed that the spin 
relaxation happens at a much faster rate than momen¬ 
tum relaxation; in case the spin-relaxation time, through 
the well-kn own D yakanov-Perel or Elliot-Yafet relaxation 
mechanism^^lllll jg comparable to the momentum relax¬ 
ation time, the overall scattering rate and mobility will 
be altered. This aspect has not been investigated here. 
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